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Abstract: Consider a large system of Brownian motions in R*^ on some fixed time 
interval [0,(3] with symmetrised initial-terminal condition. That is, for any i, the 
terminal location of the i-th motion is affixed to the initial point of the cr(i)-th motion, 
where o" is a uniformly distributed random permutation of 1, . . . , A. We integrate over 
all initial points confined in boxes A C M"^ with respect to the Lebesgue measure, and 
we divide by a normalisation (partition function). 

Such systems play an important role in quantum physics in the description of Boson 
systems at positive temperature 

In this article, we describe the large-A behaviour of the empirical path measure 
(the mean of the Dirac measures in the N paths) when A | R'^ and A/|A| p. The 
rate function is given as a variational formula involving a certain entropy functional 
and a Fenchel-Legendre transform. The entropy term governs the large-A behaviour 
of discrete shape measures of integer partitions. Any integer partition determines a 
conjugacy class of permutations of certain cycle structure. 

Depending on the dimension and the density p, there is phase transition behaviour 
for the empirical path measure. For certain parameters (high density, large time hori- 
zon) and dimensions d > 3 the empirical path measure is not supported on all paths 
[0, oo) — > R*^ which contain a bridge path of any finite multiple of the time horizon 
[0,(3]. For dimensions d= 1,2, and for small densities and small time horizon [0,(3] in 
dimensions o? > 3, the empirical path measure is supported on those paths. In the first 
regime a finite fraction of the motions lives in cycles of infinite length. 

We outline that this transition leads to an empirical path measure interpretation of 
Bose-Einstein condensation, known for systems of Bosons. 
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1. Introduction and Results 

We study the large- AT behaviour of a system of N Brownian motions B^^\ . . . , B^^^ , with time horizon 
[0,l3] in M.'^ confined in subsets Ajv C M*^, i.e., the behaviour of the system under the symmetrised 
measure 

N 

Here Sat is the set of permutations of 1, . . . , A'', jjLx]y is the Brownian bridge measure on the time 
interval [0, /3] with initial point x G Ajv and terminal point y e An and confinement to stay in A^r 
(see (2.4) below), and Z)^'""'(/3) is the normalisation 

N 

z'r\p) = 1^ E / da^i • • 7 (8) /^^^Iw )' (1-2) 

where is the set of continuous functions [0, /3] M'^. Hence, the terminal location of the i-th 
motion is affixed to the initial location of the cr(i)-th motion, where cr is a uniformly distributed 
random permutation. That is, in (1.1) we have two random mechanisms. First we pick uniformly a 
permutation and after that we pick N initial points in A^v which are permuted according to the chosen 
permutation to obtain N terminal points. Then these A'^ initial and terminal points determine the N 
random processes. Finally we average over all permutations and integrate over all initial points in the 
set Ajv- Beside the fact that the symmetrised measure Pj^'""^ is itself of interest, there are two main 
motivations in studying the symmetrised measure. 

The symmetrisation in (1.1) is described by the set of N pairs (xi, . . . ,xn; x^(^i^ , ■ ■ ■ , •'Ca-(jv)) for any 
permutation a & &n and any x G . The mixing procedure for the second entry in these pairs 
has been studied both in [DZ92] and [Tr02], which were motivated from asymptotic questions about 
exchangeable vectors of random variables. [DZ92] studies large deviations for the empirical measures 
jf X^i^i '^Yi^ where Yi, . . . , Yn have distribution Jq fi{d6) P^^ for some distribution // on some compact 
space 0, and the empirical measures are assumed to satisfy a large deviation principle under P^' for 
each 6. In [Tr02], a similar problem is studied: given a sequence of random vectors (i^/'^', . . . ,YI^^) 
such that the empirical measures J2iLi i^yev) satisfy a large deviation principle, another principle 

i 

is established for the process of empirical measures X^ll^^ ^x^^^^ where 

i 

Our second main motivation for studying the symmetrised measure Pj^'""^ stems from the applica- 
tions of Feynman-Kac formulae to express thermodynamic functions in quantum statistical mechanics. 
These thermodynamic functions are given as traces over exponentials of the Hamilton operator de- 
scribing the quantum system. There exist two kinds of elementary particles in nature, the Fermions 
and the Bosons. The state of a system of N Bosons is described by a symmetrisation procedure like in 
(1.1), whereas the state for Fermions is given with the corresponding anti-symmetrisation procedure. 
Thus one is lead to employ large deviation technique to study the large iV-limit for expectations with 
respect to the symmetrised distribution. 

Let A]sf C M"^ a sequence of subsets with A^/|AAr| ^ p G (0, oo) as iV — oo. We are going to study 
large deviations for the empirical path measure 

1 ^ 

1=1 
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which we conceive as a -measurable random measure in where Q is the set of ah continuous 

paths [0, oo) M*^ (see (2.8) below for details), under the symmetrised measure P^^™'. More precisely, 
we derive a large deviations principle for the distributions of L^v under P^^™' as A'' ^ oo (Theorem 2.2). 
(In Section 4.1 below we recall the notion of a large deviations principle.) 

To prove large deviations principles under the symmetrised measure (1.1) one cannot proceed as 
in the Gartner-Ellis-Theorem. A first approach is [AK06] and [AD06] for a different symmetrised 
measure. However, to see the peculiar correlations due to the symmetrisation one has to study the 
cycle structure for any permutation. The cycle structure allows to concatenate different Brownian 
motions to Brownian bridges with larger time horizon. Hence, a large deviations principle for the 
empirical path measure as a random element in the set of probability measures on paths with time 
horizon [0, oo) can be obtained. 

The cycle structure allows to replace the sum over permutations by a sum of integer partitions. For 
any integer A^, a partition A of A'' is the collection of integers ni > n2 > • • • > Uk > 1, k & {1, . . . , N}, 
such that Yli=i = ^- We denote the set of all partitions of N by Vn- Any partition A G Pn is 
determined by the sequence {r^}^^ of positive integers such that X^j^i kr^ = N, where we write 
rfc(A) = rfc. We call the number an occupation number of the partition. A cycle of length /c is a 
chain of permutations, such as 1 goes to 2, 2 goes to 3, 3 goes to 4, etc. until k — 1 goes to k and 
finally k goes to 1. A permutation with exactly cycles of length k is said to be of type {rk]k=i- 
Hence, each partition A G Vn corresponds to a conjugacy class of permutations, i.e., those of the same 
type, with 

N\ 

liLirklkr-^ 

elements. The conjugacy classes are the ones where the quantum mechanical trace operation is con- 
stant, see Section 4.3 for details. 

The rate function of the large deviations principle in Theorem 2.2 is given as a variational formula 
involving an entropy functional and a Fenchel-Legendre transform. The entropy term governs the 
large- AT behaviour of discrete shape measures of integer partitions. 

In Theorem 2.1 we show the thermodynamic limit of the normalisation Z^'^\p), and provide a 
variational formula for it. We analyse the variational problem, i.e., the minimisers for the entropy 
functional, in Theorem 2.5. In particular, we derive a variational expression for the specific free energy. 
A phase transition, i.e., a singularity of the free energy, exists for dimensions d >3 depending on the 
density p and f3. For high density, respectively for long time horizons [0,/?], the specific free energy 
is independent of the density. This is the Bose-Einstein condensation for a non-interacting system of 
Bosons in the thermodynamic limit (see Appendix 4.2). 

Our main results, the large deviations principle for L^r under Pj^^™', in Theorem 2.2 and the analysis 
of the rate function in Theorem 2.6, give a phase transition for the empirical path measure. Let the 
set 

where is the set of paths [0, and (8)fe/3 is the splice of two paths, be given. The empirical 

path measure for dimensions d = 1,2, and p < oo, or p < pc for d> 3 has support on paths in any 
Af;, where one can insert any finite number of Brownian motions with time horizon [0, /3], i.e., for any 
A; G N one can concatenate exactly k Brownian motions to paths u e A^. This is due to the cycle 
structure of the permutations and the Lebesgue integration over all initial points in the definition of 
the symmetrised measure P^^™'. If the density p is high enough for d > 3, i.e., p > pc (or equivalently, 
if the time horizon is sufficiently large for given density, i.e., /3 > /3c, for d > 3), the mean path measure 
has positive weight on paths with an infinite time horizon, that is, concatenation of any finite number 
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of Brownian motions with time horizon [0, /9] to obtain paths in the A^s, is not sufficient. Then an 
excess density [p — Pc) of Brownian motions with time horizon [0, [3] concatenate to infinite long cycles. 

Hence, we have an empirical path measure interpretation of Bose-Einstein condensation. This novel 
interpretation adds to the existing notions of Bose-Einstein condensation ([FW71],[Sii02],[Uel06]), 
and allows to study systems of interacting Brownian motions. Future work will be devoted to this 
case [A07a]. Interacting Brownian motions in trap potentials have been so far analysed without 
symmetrisation, in particular, finite systems for vanishing temperature in [ABKOGa] and large systems 
of interacting motions for fixed positive temperature in [ABK06b]. 

Let us make further remarks on related literature. 

An important work combining combinatorics and large deviations for symmetrised measures is 
[T6th90]. Toth [T6th90] considers N continuous-time simple random walks on a complete graph 

with pN vertices, where p G (0, 1) is fixed. He looks at the symmetrised distribution as in (1.1) 
and adds an exclusion constraint: there is no collision of any two particles during the time interval 
[0, (3\ . The combinatorial structure of this model enabled him to express the free energy in terms 
of a cleverly chosen Markov process on Nq. Using Preidlin-Wentzell theory, he derives an explicit 
formula for the large-iV asymptotic of the free energy; in particular he obtains a phase-transition, 
called Bose-Einstein-condensation, for large /3 and sufficiently large p. 

Scaling limits for shape measures of integer partitions in Vn under uniform distribution are obtained 
in [Ver96] . Large deviations from this limit behaviour are in [DVZOO] , where large deviations principles 
for scaled shape measures for partitions as well for strict partitions under uniform distributions are 
derived. Motivated by the statistics of combinatorial partitions, illustrated by Vershik in [Ver96], 
Benfatto et al. derived in [BCMP05] limit theorems for statistics of combinatorial partitions for the 
case of a mean field Bose gas in the grandcanonical ensemble. Here, in contrast to the canonical 
ensemble, only the mean of the particle number is fixed. [BCMP05] uses Fourier analysis of the 
corresponding traces to derive a complete description of the statistics of short and long cycles. For a 
perturbed mean-field model the density of long cycles for a perturbed mean-field model is analysed in 
[DMP05]. 

In Section 2 we describe our main results, where we comment on the empirical path measure 
interpretation of Bose-Einstein condensation at the end of that section. Section 3 is devoted to the 
proofs, and in the appendices in Section 4 we provide the notion of large deviations principle, and 
review facts on the zeta function of Riemann and on quantum statistical mechanics for systems of 
Bosons. 

2. Results 

In this section we present our results. The phase transition result in Theorem 2.6 leads us to a 
probabilistic interpretation of Bose-Einstein condensation on the path measure level. I Throughout 
the paper, we fix /3 > 0. Let a:={uj: [0, oo) ^W^-.u) continuous} be the set of continuous functions 
[0, oo) — W^. The elements in O are called trajectories or paths and we denote by = {w : [0, kf3] — > 
W^: uj continuous}, /c G N, the set of paths for time horizon [0, k(5]. We write Q,p for We equip J7 
(respectively Jlfe) with the topology of uniform convergence and with the corresponding Borel cr-field 
B (respectively B^). We consider Brownian motions, B'--^\ . . . ,B'-'^\ with time horizon [0,/?] as 
random variables taking values in For the reader's convenience, we repeat the definition of a 
Brownian bridge measure; see the Appendix in [Sz98]. We decided to work with Brownian motions 
having generator A instead of ^A. We write P^^ for the probability measure under which B = B^^^ 
starts from x G M"^. The canonical (non-normalised) Brownian bridge measure on the time interval 
[0, P] with initial site x G M'^ and terminal site y G R'^ is defined as Hx,yiA) = Pj;(S G A; Bp G dy)/dy 
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for A C measurable. Hence, the Brownian bridge measure for a Brownian bridge confined to a 
subset Ajv C M*^ is defined by 

P'xlyi^) = ) AcClp measurable. (2.4) 

Then i2x]y is a regular Borel measure on O^. The normalised Brownian bridge measure is defined as 

Pjf = <f/<f(f^/3), (2.5) 

which is a probability measure on Qp. 

We conceive the empirical path measure as a random element in A^i(O), hence, we need a convenient 
extension of any continuous path [0,(3] — M'' to a continuous path [0, oo) — ^ in the definition (2.8) 
of the empirical path measure. For any x G we denote by the Brownian probability measure on 
O, i.e., the canonical Wiener measure with deterministic start in a; G M'^ ([Gin70]). In the following 
we write alternatively ut or u;{t) for any point of a path u. Given a path a; G with time horizon 
[0, P] define 

= Su. ®p P^'^"^ e Mi{Q,B), (2.6) 

where the product (X"^ is defined for the "splice" of two paths, i.e., for lo G and G 17 define aJ G O 
by a;(t) = (jo{t A /3), t G [0, oo), and u <Si/3 w G such that uj (g)^ a3 = tD if u;(0) 7^ '^(/3) and 

r,m-/ fortG[0,/3] 
a;®;ja;(t)-| for t G (/3, 00) ^^'^^ 

if a;(0) = u}{P). The mapping w G i7/3 1-^ P^,''' G A1i(r2,;S) is measurable, and the family {P^,''^: a; G 
Qfj} satisfies the Markov property ([DSOl, Lemma 4.4.21]). Hence, the empirical path measure 

1 ^ (■) 

i=l 

is ^2^^ measurable. Here co = (w'^V • • ) '^*^') G Our main result concerns a large deviations 



principle for the distributions of Ljv under the symmetrised measure . Recall that P]^^™ is a 



probability measure on 

To formulate the rate functions we need some notations. 
Now we introduce the rate function. Let 

M = {Qg [0,lf- ^Q{1) = 1,Q(0 > Qil + 1)V/ G N} 

be the set of monotonously non-increasing probability functions on N. For Q e M define Q{k) = 
Q{k) - Q{k + 1) for any A; G N. For d > 1 let 

be given, and define the functional 

S{Q)=f^Q{k)[log-S^-l), QGM. (2.10) 
fe=i ^ ^'^) 

The rate function is given by 

/(--)()") = Jnf^ {S{Q) + I^^Kfi)} - X(/?, p), fi G Mi{n), (2.11) 

where 
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sup |(F,M)-^Q(fc)logEfo(e^(^))|, i^eMiin), (2.12) 

and where the function x{P:P) •= infQex{«5(Q)} is given as the negative logarithmic Umit of the 
partition function see Theorem 2.1, and where Cb(fi) is the space of continuous bounded 

functions of the paths in Q. Eq q denotes the expectation with respect to the Brownian bridge proba- 

biUty measure Pq'^ defined in (2.5) extended as in (2.6) to a probability measure in A4i{i}). Here, Z^*^^ 
is a Fenchel-Legendre transform, but not the one of a logarithmic moment generating function of any 
random variable. In particular, and therefore also /(^y™', are nonnegative, and I'-'^^ is convex as a 
supremum of linear functions. There seems to be no way to represent /("^'(/x) as the relative entropy 
of n with respect to any measure. 

We need the following limit for the partition function to derive the large deviations principle. 

Theorem 2.1 (Limit of the partition function). Let p G (0, oo) and An C with An T o-f^d 
A'"/|Ajv| —>■ p as N ^ oo and with dist(0, OAn) = cN for some c > 0. Then 

This Theorem gives the normalisation part xif^^p) = ^^^QeMi'^iQ)} the definition of the rate 
function in (2.11) for the large deviations principle in the following Theorem 2.2. 

Theorem 2.2 (Large deviations for Ln). Let p G (0, cx)) and An C with An T ffi*^ and 
N/\An\ p as N ^ oo and with dist(0, OAn) = cN for some c > 0. 

Under the symmetrised measure P^^™' the empirical path measures {LN)Nen satisfy a large devia- 
tions principle on Mi{^) with speed N and rate function 

The proof of Theorem 2.2 is in Subsection 3.3. The proof is based on the following variational 
formula for the logarithmic limit of the normalisation in Theorem 2.1 and on the analysis of the rate 
function in Theorem 2.5 below. 

We give a brief informal interpretation of the shape of the rate functions in (2.11) and (2.12). As 
remarked earlier, the symmetrised measure P^^™^ arises from a two-step probability mechanism. This 
is reflected in the representation of the rate function in (2.11): in a peculiar way (we describe it 

at the beginning of Section 3), the term S{Q) — xiPiP) describes the large deviations of the discrete 
empirical shape measure for integer partitions. The discrete empirical shape measure Qn, defined as 

1 ^ 

QN:PN^Miif^),X^-Y.rk{\), (2.14) 

k=- 

governs a particular distribution of A'" independent, but not identically distributed, Brownian bridges. 
Under this distribution, Ln satisfies a large deviations principle with rate function which can also 
be guessed from the Gartner-Ellis theorem [DZ98, Th. 4.5.20]. The presence of a two-step mechanism 
makes it impossible to apply this theorem directly to P^^™'. 

Let us contrast this to the case of i.i.d. Brownian bridges B'-^\ . . . , B'-^\ starting in the origin, i.e., 
we replace Pj^^™' by (Pq q)®^. Here the empirical path measure Ln satisfies a large deviations principle 
with rate function 

/(m)= sup {(F,M)-logEjo(e^(^))|, 

as follows from an application of Cramer's theorem [DZ98, Theorem 6.1.3]. Note that L(p) is the 
relative entropy of p with respect to Pq q o B~^. Although there is apparently no reason to expect a 
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direct comparison between the distributions of Ltv under P*^^""^ and under (Pq q)^^, the rate functions 
admit a simple relation: it is easy to see that I'-'^^ > I for the measure Q e M with Q{k) = di, since 

oo 

-^Q(A;)logESJ[e^(^)] > - log E (e^(«) ) . (2.15) 
fe=i 

In particular, I^'^'"") > /. 

Remark 2.3. The techniques of the proof of Theorem 2.2 apply also to a proof of a large deviations 
principle under the symmetrised measure P^^™' for the empirical path measure Ljv = ^/^^f=\ <^b(<); 
which is a random element in M.i{ytp). The rate function is 

oo 

= inf {S{Q) - sup - ^ Q(fc) logE^J (e^'=° 1 1 e >fi(%). 

Similar results hold for the mean of the occupation Yj^ measures, 



However, these rate functions seem not to give enough information to derive the phase transition as 
in Theorem 2.6, and to obtain a probabilistic interpretation of Bose-Einstein condensation. 

Our large deviations result is accompanied by an analysis of the variational formula for the rate 
function (2.11), i.e., the analysis for zeros of the rate function. This gives the proof phase 

transition for empirical path measures depending on the dimension and the density parameter in 
Theorem 2.6. We analyse first the variational formula for the limit of the partition function. 

Remark 2.4 (Free energy). The variational formula for the limit in (2.13) gives the specific free 

energy f(P,p) := limjv-»oo ~l//5|-^Ar|-Z^^''™^(/3) for inverse temperature (3 and density p, i.e., 

This is the specific free energy for an infinite systems of non-interacting Bosons in the thermodynamic 
limit, i.e., the limit where N/\A.n\ — ^ p as N ^ oo. 

We analyse the variational formula for %(/?, p) (2.13), and we derive an expression for the free energy 
/ as a function of (3 and p. Define a dimension dependent critical density 



d)^C(f), ford>3 



d\ s-^ , _d 
k 2. 



p^=.l (47r/3)'*/2H2; ' - , (2.17) 

[ +00 , for d = 1, 2 

where ( is the zeta function of Riemann, 

k=l 

Furthermore, denote by gs{ct) the so-called Bose functions (see (4.57) in Appendix 4.2) 

oo 

gg[a) = ^ k~^e~°''' for aU a > and ah s > 0. 
k=l 

For any p < pc we denote by a = «(/?, p) the unique root of 

^ oo 

p = — y fe-'^/2e-°^ (2.18) 
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The essential difference in d > 3 and d = 1,2 lies in the fact that in the latter two cases the corre- 
sponding Bose functions, gi{a) respectively gi{a), diverge as a ^ (see Appendix 4.2 and [Gra25]). 
For d= 1,2 there is a unique a for any density p < oo. For d>3 there is such an unique a given only 
for densities p < pc- Hence, this is the mathematical origin of the so-called Bose-condensation (see 
Appendix 4.3), where for d > 3 and p > pc particles condense in the zero mode state. This corresponds 
to the phase transition for the empirical path measure, when the empirical path measure has positive 
mass on paths which do not return to the origin at any times which are an integer multiple of /3. 

In the following theorem we analyse the variational problem (2.13). 

Theorem 2.5 (Analysis of the variational formula for xiPiP))- ^or any p < oo in dimensions 
d = 1,2, and p < pc in dimensions d > S, there is a unique minimiser Q E M of the variational 
formula (2.13) with 

Qik) = ^ \2 for ken, (2.19) 

p(4^/3)°'/2fci+2 

whereas for dimensions d > 3 and densities p > pc, there is no minimiser for the variational problem 
(2.13), but the infimum is attained for any minimising sequence {Qn)n£N of Qn £ M such that 
Qn Q* as n ^ OO. 

The specific free energy for d> S is given by 

~ i4nf3r/2f3 9d±?.{a) - ^poi , for p< Pc 



f{p,p) 

and for d = 1,2 by 



(2.20) 



where a is the unique root of (2.18). 

This result leads to the analysis of the rate function . Let 

= {a;(8)fc^C: a; G J^fc,a;(0) =uj{k^),$, en} cn,k eN, (2.22) 

be the set of paths in J7 which result from the splice (2.7) of Brownian bridges paths of time horizon 
[0, kp] with any path ^ G Jl. 

Theorem 2.6 (Analysis of the rate function /('^y™)). Under the assumptions of Theorem 2.2 the 
following holds: 

(i) d = 1,2. A unique minimiser p,* G Mi{n) of the rate function p i— /('''""'(/x) is given with 

(ii) d > 3 and p < pc- A unique minimiser p* G of the rate function p i-^ /'"^""'(/x) is given 
with Y^ken kp*iAk) = 1. 

For p > Pc there is no unique minimiser be given, but there exist minimising sequences 
iPn)n>i, pn £ ■Mii^,), with XlJ^i ^/"n(^fc) = 1 for any n G N such that pn ^ fJ-^ ^ Mi{fl) 
weakly as n ^ oo with X^^i kp^{Ak) < 1. 

Proof. := S{Q) - xiP, p) > for all Q e M. Moreover, I'«)(/x) > for any Q G X, which can 
be seen putting F = in the supremum (2.12). The functional <S can be written as 

S{Q)=qH{P\P*) + qlog-^-q ,Q E M, 

where q = Ylk^=i Q{k)i 1* = YlV=i ^■^d H is the relative entropy of the two probability measures 

P = q~^Q and P* = {q*)~^Q*. Hence, the level sets oi Q i-^ ^xiQ) compact. 
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(i) Case d = 1,2. Assume that fi is a zero of /'^y™'. Since the level sets are compact, there is aQ^ £ M 
that minimises the formula on the right hand side in the definition of /(''y™' (2.11). Clearly /C^'') < oo. 
As both parts of the rate function are positive, we conclude that Q'^ is the unique minimiser of 

the functional S (see Subsection 3.2 in the proof of Theorem 2.5), i.e., S-)^{Q'^) = 0. Note further that 
a subclass of Cb(^^) consists of those function Fi;,k eN, such that 



rkl3 

Fk{i^ (^kp 0=1 /('^s)ds for ah ^ G ri, 

^0 



for continuous bounded functions fk'^'^ R The Euler-Lagrange equations yield the following 
system of equations 



where H = Fk,k e N, is any of the functions defined above. F = solves this system, and if one 
chooses /fc = 1//?, A; G N, in the definition of the functional Fj^, one obtains that 

kn{Ak) = kQ^{k) ,keN. 

This identifies fi as /i*, since YlT=i kQ^{k) = YlT=i Qi.^) = 1- 

(ii) Case d > 3. For densities p < Pc the assertion follows as in (i). Let p > pc- Define the sequence 

{Qn)n&N,Qn & M, by 

where Q* is given in (2.9) (see Subsection 3.1), i.e., 

p(47r/3)'^/2A;i+'i/2 

Clearly, Qn Q* strongly as n ^ oo, and lim.n^^ S^{Qn) = 0. As 

oo oo ^ ^ 

the infimum is not attained in M. Hence, from the considerations in (i), we get a sequence {pn)n£N of 
zeros iJ,n G M-i{Q) of the rate function JC^") with /x„(^fe) = Qn{k) and Yl'k^=i kpni^k) = 1 such that 
/('3n)(^^) = 0. Then, 

/'=^"^'(Mn) < inf {s^iQ) + I^^'if^n)} < S^iQn), 
QeM I- J 

and from the upper semi continuity of the functional 5, the strong convergence of Q„ ^ Q* as n ^ oo, 
and from the fact that the level sets of 1^"^'^'^ are compact, one derives that pn —>■ P'^ & Mi{^) weakly 
as n — *■ oo, and lim„^oo -^*°^™'(/"n) = 0- ^ 
Let us draw an easy corollary from this theorem. 

Corollary 2.7 (Law of large numbers). Under the assumptions of Theorem 2.6 the following holds. 



(i) For d = 1,2, and any density p < oo, there is a law of large numbers. Under the probability 
measure the sequence (Lat )jvgn converges in distribution to the measure p* G ^Al{Q). 

ii) For d > 3 and p < Pc there is a law of large numbers. Under the probability measure P^^™', 
the sequence (Ljv)iveN converges in distribution to the measure p* G A^i(O). 
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The main conclusion of the large deviations principle in Theorem 2.2 and Theorem 2.6 is the follow- 
ing phase transition for the mean empirical path measure, which gives a path measure interpretation 
of Bose-Einstein condensation (BEC). 

Path measures and their interpretation as Bose-Einstein condensation 

Let N Brownian motions with time horizon [0, (3] confined in sets Ajv C M*^ given such that Ajv T 1^*^ 
and dist(0,5Ajv) = cN for some c > and A^/|AAr| p E (0, oo) as AT — oo. Then the following 
holds: 

(i) For P > there is a, pc = PciPi d) such that: 

no BEC: Case p < pc for d > 3, p > fov d = 1, 2: 

Ln^ p* € Mi{n) under P^^""^ as AT ^ oo with kp*{Ak) = 1 

BEC: Case p < Pc and d> 3: 

Ln ^ P^ e Mi{n) under P^'""' as AT ^ oo with YlT=i kp^i^k) < 1- 

(ii) For p G (0, oo) there exists a 



such that: 

no BEC: Case /3 < /3c for d > 3 and /? > for d = 1, 2: 

Ln p* e Mi{n) under P^^""' as iV ^ oo with kp*{Ak) = 1 

BEC: Case /3 > /3c and d > 3: 

Ln^ p^ e Mi{d) under P^^"") as AT ^ oo with ^^^^ hp^{Ak) < 1. 

If d = 1, 2, or p < Pc for d> 3, the mean empirical path measure has support on those paths in which 
one can insert, starting from time origin, a concatenation of any finite number of Brownian motions 
with time horizon [0, f3], i.e., for any A; G N one can find in paths uj^ G Ak exactly k Brownian motions 
concatenated to a Brownian bridge with horizon [0,/c/3]. This follows from the concatenation of the 
Brownian motions due to the cycle structure of the permutations and due to the Lebesgue integration 
of any initial position in the definition of the symmetrised measure P^^""^ (1.1). If the density p is high 
enough for d > 3, i.e., p > pc (or equivalently, if the inverse temperature is sufficiently large for given 
density, i.e., /3 > /?c, for d > 3), the mean path measure has positive weight for paths with an infinite 
time horizon, that is, concatenation of any finite number of Brownian motions with time horizon [0, f3], 
i.e., any finite cycle path in A/^., is not sufficient, because there is an excess density {p — Pc) of Brownian 
motions with time horizon [0, (3]. These motions concatenate to infinite long cycle, that is, these cycles 
grow with the system size in the thermodynamic limit. The fraction of these motions is 



Let us briefly contrast our results to the ones derived in [AD06] and [AK06]. There the empirical 
path measure was studied under the symmetrised measure 



where m G A^i(M°') is an initial distribution for the Brownian motions. Multilevel large deviations 
principle are obtained, however, a phase transition for the mean path measure does not exist. Only 
in the special case where m is the Lebesgue measure restricted to a fixed box A C M*^, it turns out 
that the rate function for finite time horizon in the limit A^ ^ oo is precisely the Donsker-Varadhan 
rate function. This rate function appears in the limit of infinite time horizon, hence its appearance is 
an indication of the presence of infinite long cycles. Note, this corresponds to the limit A^/|A| oo. 





(2.24) 
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i.e., infinite density. In [AD06] a similar result was shown for quantum spin system with mean-field 

interaction. 

Our results here arc essentially different, because of the definition of the measure P^^™' and of the 
limit N/\A]\[\ ^ p G (0,oo), and go far beyond the ones in [AD06] and [AK06]. 

Finally, let us remark that our results show that the discrete empirical shape measure 

1 ^ 
fe=- 

shall satisfies a large deviations principle under the distributions /xjv given by 

The distribution /Xjv is not the uniform distribution on the partitions Vn (d > 3). Limit theorems 
for scaled shape measures have been obtained in [Ver96] for uniform distributions of partitions. Large 
deviations under the uniform distribution from this limit behaviour with appropriate scaling of the 
shape measures have been derived in [DVZOO]. In [A07b] we will study such large deviations for the 
non-uniform distributions given above. 



3. Proofs 



In this section we proof our the Theorems 2.1, 2.5 and 2.2. In Subsection 3.1 we prove the limit of the 
partition function, i.e., the thermodynamic limit of the free energy and its expression through a varia- 
tional formula for discrete shape measures. The large deviations principle is proved in Subsection 3.3, 
and the analysis of the variational formula is in Subsection 3.2. 



3.1. Proof of Theorem 2.1. The key idea is to use the cycle representation for permutations and 
to replace the sum over permutations in the definition (1.2) of the partition function by a sum over 
integer partitions of N. For any partition A G Vn 

we assign the discrete empirical shape measure 
for the occupation numbers {r^}^^ with X^^i krk = N, i.e., the probability measure G A^i(N), 
defined as 

Qn:Vn^Mi{N) 



^ (3.25) 



A-Q^(-) = ^E^fc(A). 

fe=- 

Q^{1) = for / > AT. Here Q^(l) = J2k= 

^ Tfc is the number of components/ cycles. The mapping 
Qjv is injective, which can easily be seen as follows. Qp^ = for A, A' G Vn with A 7^ A' implies 
successively r7v(A) = rjv(A'), . . . ri(A) = ri(A') and hence A = A' in contradiction with our assumption. 
For any partition A G Vn, i-e., for any set of occupation numbers {rk}^^i with "^^^i kr^ = N, we 
have A^!/ rifcLi permutations of type {r^j^j^. For a given permutation A G Vn we regroup 

the product of the Brownian bridges as follows. In the following we write for rfe(A). For any 
A; G {1, . . . , N} with > 1 we have exactly groups/cycles of length k which involve altogether kvk 
Brownian motions, 

( dxi • • • dxk ® /x£3 ® • • • /xff,^) ( j^^ dxk+, ■ ■ ■ dx,k Mfi,.,^, ® • • • ® mK,) 
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numbered here from 1 to krk- The stxTicture of the product measure in each of the groups/cycles 
is identical, hence we define the measure 



(3.26) 

kf3,N 



where we used spatial shift invariance of the Brownian motion and concatenated the Brownian bridges 
to a Brownian bridge of time length kp starting and terminating at the origin. Here, (see (2.4)), 

^.j.^^) ^ MA,B,,e dy,B,,,,,,cA.)) ^^^^^ ^e,,,,,ble. 

This concatenation is possible because of the indicator to stay within Ajv- Hence, we obtain for the 
partition function 

AT 



Z'r' = 1^ E / dxi • • • / da:^^ (g) {^^ ) 

En(Mr)(®(''Srr*)(%"). 



AG-Pat fc=l ^ '^^ k=l 

where we replaced the sum over permutations by a sum over partitions. To study the large A'^- 
behaviour of the partition function we rewrite the expression on the right hand side of the latter 
expression. We sum over all elements in the image 

N 

Mn := Qn{Vn) =C {Q G {0, l/AT, . . . , 1}^ : ^ Q{1) = 1, Q{1) > Q{1 + 1) VZ e N} 

1=1 

of the discrete empirical shape measure Qn, which are probability measures Q: N ^ {0, 1/iV, . . . , 1} 
with Q{1) > Q{1 + 1) for all / G N and X;/=i <9(0 = 1 and Q{1) = for all / > N. The occupation 
numbers for a partition A G Vn are then given by = N{Qj^{k) — QN{k + 1)), k = 1, . . . , N , where 
Qn = Qn- We write Qwik) = QN{k) — QN{k + 1) ioi k = 1, N, in the following. Thus we obtain 

On the right hand side we have to evaluate the probability mass on . Now we use that /j^q^q^ , k G 
{!,..., N}, is a regular Borel measure on flk with total mass equal to 

MoJ'''(^^fe)=Pfc/3,A.(0,0). (3.27) 
The transition density Pk/3,AM is defined as 

Pkp,ANi^^y) = -^^^^^^^^—^ -fJ^NikP,x,y), x,yeAN, (3.28) 

with 

rAN{k(3,x,y) = E^{taj, < k^; (Br^N^y)) (3.29) 

and 

Pkp{x,y) = ^^^^'^^^ = {4nPk)-<^/'e-\^-y\y^^^ 
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where ta^ is the exit time of the Brownian motion of the set An ([CZ95]). Without loss of generahty we 
let An C M.'^ be centered such that dist(0, SAjv) = aN for some a > 0. For any c with < kp < c^/d 
and \x — y\ < c< dist(x, ^A^v) A dist(y, OAn) one gets ([CZ95]) the estimation 

rA^{kp,x,y) < (47r/3fc)-'^/2e-^. 
Applying this with c = \/d(3N we get 

(47r/3fc)-''/2(l - e-<^W) < fj'^^f {n^) < {ATrPk^/^ 

Hence, 

^ I A I \'0(k) 
(1 _ ^-dN/Ap^N y / TT _^ _\ ^ z'^yn.)^^^ 



Recah from (2.10) that 



~ „^ VJ--L (Nn(h\\\hNQ(k)(A^RU\d/2NQ(k)) 



^ - Q(A;) 



5(g) = ^Q(fe) logf^-1 



fe=i 



for Q e Mn- 



Stirling's formula N\ Ri N^+^/'^e-^V2n{l + 1/127V + 0{N^)) gives the following upper and lower 
bound 



/ ^ ^NQ{k) , 

< C-^-gATioga+oCAT)) ^ exp(-iV5(g)), 

and 

Z%'-\P) > (1 - e-'^W)JVe-iiv« iogiVg-iviog(i+0(iV)) ^ ^^^g)^ ^ ^^^^^^ 

for a G (^)l)- The error bounds in (3.30) and (3.31) due to Stirling's formula follow from the 
observation that there is a a G (^, 1) such 

AT 



: NQ{k) > 1, ^ kr^ = N} = N'^. 



k=l 

Recall that NQ{k) = rk,k = 1,...,N, the existence of such an a can be seen from 

N N 

'Aik : NQ{k) > 1, fcrfc = iV} < : NQ{k) = l,^krk = N} < cN^ . 

k=l k=l 

Thus an upper (3.30) and a lower bound (3.31) is established. To proof Theorem 2.1 we need to 
study the large A^-limit for the upper and lower bound, and the continuity of the functional S. 

Proposition 3.1. Let (Ajv)jv>i be a sequence of centred cubes with An | and A''/|Ajv| —>■ p as 
N ^oo. 
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(a) Assume additionally that dist(0, 5Ajv) = cN for some c > 0. Then 

liminfllog J2 exp(-iV5(Q)) >- inf 5(Q). 
Q&Mn 

(b) 

lim sup 1 log V exp ( - NS{Q)) < - inf S{Q). 

Proof. (a) It suffices to construct, for any Q G Al, some sequence {Qn)n>i with Qn e Mn which 
converges strongly to Q as A'" — oo such that 

limsup5(Qjv) < S{Q). (3.32) 

AT— >oo 

Let Q e M. with J2k£n Qi^) = 1 be given. Recall that [xj is the largest integer smaller or equal to 
xeR. For a G (0, 1) define 

i^^, ifA; = 2,...,LiV-J ■ (3.33) 
, otherwise 



Qjv(A;) = < 



It is clear that X;f=i QN{k) = 1, and Qn{1) > Qn{1 + 1) for i > 2 and g7v(l) > 1 - Efe>2<3(^) ^ 
Q(l) > Q(2) > QnC^) follows from the fact that x <y implies [x\ < [y\ for x,y e M+. The sequence 
(QAr)AreN converges strongly against Q as AT — oo because 

^\Q{k)-QN{k)\=Y.\Q{k)-QN{k)\+ Yl Q{k)<N''-'+ Yl ^(^)- 

k=l k=l A:>[7V"J A;>[Af«J 

To prove (3.32) note that the functional S is lower semi continuous, because S is the sum of the relative 
entropy of Qjv with respect to Q* restricted to {!,..., [A^^J} and the total mass Qwik) < 1, 

and the relative entropy is lower semi continuous. Hence, liminfjv-»oo«5((5Ar) > S{Q), because of the 
strong convergence of Qtv to Q as A'^ ^ oo. We are left to prove the upper semi continuity (3.32) of 
the functional 5 for the sequence {Qn)n>i- This is seen from 

TV 

s{Qn) = Yl {Q^(^) log ^^(^) - QN{k){i + log g*(fc))) , 

k=l 

where we have to analyse only the entropy term, i.e., the first term on the right hand side. The entropy 
is defined by H{Q) = - Q{k) log Q{k) for Q e M. 

H{Qm)-H{Q)= ^Q^(A;)log^^+ ^(gjv(A;)-Q(A;))logQ(A;) 

k=i ^ QW k=i ^334^ 

+ Y Q{k)iogQ{k) 
fc=[Ar«j+i 

In the first two terms on the right hand side of (3.34) it is sufficient to consider only sum indices 

k £ {1, . . . , [A^"J} with Q{k) > \/N . The terms for /c = 1 in the two sums on the right hand side 
vanish due to the strong convergence. The remaining first term on the right hand side of (3.34) is 
bounded by 

[AT" J 

Y \QN{k) - Q{k)\ < 2Ar"-i ^ as AT ^ oo, 

k=2 
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the second term is bounded 



J2 {Q{k)N-Q{k))\ogQik) 



fc=i 

Q(k)>l/N 



< 2iV"-^ logiV ^ as AT ^ oo. 



The third term clearly vanishes in the limit N ^ oo. The latter estimations together with (3.34) 
finish the proof for the upper semi continuity, i.e., limswpj^^^ S{Qn) < S{Q), and therefore (3.32) 
is proven. This together with the lower bound in (3.31), i.e., all the error terms vanish in the limit, 

proves (a). 

(b) For the upper bound we need to estimate the cardinality of the set Mn- The set Mn consists of 
all antitone mappings Q: {1,...,N} ^ {0,1/N,2/N,...,1} with E£iQ(0 = 1 which are elements 
in the image Qn{Vn) of the discrete empirical shape measure Qn (3.25). For A G Vn there are at 
most A^2 indices k such that the occupation number for a component/cycle of length k is nonzero, i.e., 
''fc > 1- These indices are called occupied. One maximises the number of occupied indices if one starts 
with index 1 and fills each component/cycle only ones. Then, at index [A^^/^] we already exceed the 
total number N, because X^fcLi kr^ = N. Here \x\ is the smallest integer greater or equal to x G M+. 
Hence, there is a a G (^, 1) such that 

^{k G {1, . . . , A^} : rfc(A) > 1} = tl{A: G {1, . . . , A^} : NQ^{k) > 1} < AT" for A G Vn- 

We choose \N°''\ indices out from N and we denote this set by M. M has cardinality [AT"] . The 
cardinality of the set {/ : M ^ {l/N, 2/N, ...,!}: /(a) > f{b) for a <b,a,b e M} is then given by 

: M ^ {1/N, 2/N, ...,!}: /(a) > f{b) for a < 6, a, 5 G m} = + 1) ■ ■ .(N+ \N-] - 1) 

_{N+ [A^"] - 1)! 
~ (A^- 1)! [Ar«]! ■ 

(3.35) 

Recall that any clement in the image Qn{T^n) is antitone and a probability measure, whereas the 
latter condition is relaxed in (3.35). We prove (3.35). Let the set M = {mi, . . . ,m|-7V"]} be given with 
the natural order. Then the mapping 6162 • ■■^n'-^ + - 1), ^'2 + ([A^"] -2),. . . ,6pArc]} 

is a bijection between the set of all antitone mappings f : M ^ {1/-^! 2/A/", . . . , 1} and the set of all 
I" TV"] -subsets of {1, 2, . . . , [A^"] + A^ — 1}. Here wc used the word representation of the mapping $. 
A^(A'^ + !)••• (AT + [A'^"] — 1) counts the number of sortings of [A^"] objects into A'' linearly ordered 
boxes. 

The occupation numbers = for all the remaining indices /c G {1, . . . , A''} \ M. Thus Q%{1) = 
for all Z G {1, . . . , A/"} \ M with I > maxM. Otherwise Q^{1) is constant determined by the indices 
of the set M, i.e., Q^{1) = l/-^Z^feGM,fc>i''feW \ M and I < maxM. Hence, to 

estimate the cardinality of the image set A4^i wc have (|-jyL^) ways to choose a subset M C {1, . . . , A^} 
of cardinality jjM = [A^"] . The remaining choices for the image elements are then given by (3.35) 
because all the indices in the complement {1, . . . , A'^} \ M are determined by the indices in M. Hence, 
using Stirling's formula again, we get 

" ^ - \\m]J (AT- 1)! [AT"]! - \ N J ^ ^ 

for some C > 0. 

Let Qn G Mn be minimal for S on the set Mn and let Q G denote the limit limjv->oo Qiv- 
From the upper bound (3.30) we obtain 

Z(^-'(/?) < CA^2riV«l(l-a)(^i _ ^^(g^)^_ (3.37) 
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By the continuity of Q i— > S{Q), i.e., by the lower semi continuity of S, the assertion (b) of Proposi- 
tion 3.1 follows. 

□ 



3.2. Proof of Theorem 2.5. We analyse the variational formula in (2.13) for the functional S and 
prove Theorem 2.5. Note that any montonuously non increasing probability measure Q € is in 
one to one correspondence to Q G {Q : N ^ [0, 1] : Z^^Li Q(^) < !}• If Q,P ^ M the equality 
Q = P implies Q{1) - Q{n) = P{1) - P{n) for n G N, which implies Q{1) = P{1) and hence P = Q. 
Furthermore note that the functional S is convex and lower semi continuous and can be written as 

S{Q) = qH{P\P*) + g log ^ - g, (3.38) 

where q = J2T=i Q{^)^Q* = Ylk^i Q*{k) and H the relative entropy of the two probability measures 
P = <r^Q and P* = {q*)''^Q*. A glance at the first term in (3.38) would imply that P = P*, because 
the relative entropy is zero if and only if P = P*. 

For d= 1,2, and any density p < oo, there is a unique solution a of (2.18). Hence, Q e M with 

^ g-afe 

Q(k) = J for ken, 

is the unique solution of the variational formula (2.13) with Yl'k^=i Q{k) = 1- The free energy is given 
by p/f3S{Q), and (2.21) is proved. 

We turn to the case (i > 3 in the following. We fix p < p^. Minimising only the relative entropy in 
(3.38) is not sufficient because Q = Q* would imply that 



This follows from the fixed density in (2.18), i.e., p = ^^^^■^a/2 9d/2ioi) and 5^/2(0;) < 5d/2(0) (see e.g. 
[Gra25]). Hence the unique solution for the variational problem in (2.13) is given hy Q e M with 

„— afe 

Q(k) = J for keN, (3.39) 

p(47r/?)'^/2fci+i 

and the free energy is given as 

^ E fci+4(4,^)./2(-"^ - 1) = -^(i^^^ --fP 

for p < pc- 

For the case p > Pc define the sequence ((5n)neN by 
Recall that Q„(0 = Q„(l) - EL=\Qn(^) for 

any I G N. Then Qn G M.,n G N, because Qn{k) > 

Qn{k + l),ke N, and 



fc=l fc = l 
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We have Qn Q* strongly as n — > oo and for any n G N the evaluation for the functional is given 
by 



s(Q„) = sm - + (g-(„) + <i^) log (^^t^ 

np \ np J \ Q*{n) 



np 



Hence, lim„^oo S{Qn) = S{Q*). The infimum is therefore not attained within the set for p > pc, 
only minimising sequences do exist. 

3.3. Proof of Theorem 2.2. The following is a reformulation of Theorem 2.2, making explicit what 
a large deviations principle is. Let Aat C with dist(0, ^Ajv) = cN for some c > and A'"/|Ajv| — > 
p G (0, oo) and Ajv T l^'' as AT — oo, then Theorem 2.2 is equivalent to (i)-(iii): 

(i) Lower bound of the large deviations principle 

For any open set G C A^i(r2), 

1 



liminf ^ logP^^'"'(L7v G G) > - inf f^^'^^ii). (3.41) 

(ii) Upper bound of the large deviations principle 

For any compact set K C Ali(ri), 

lim sup 4 log F^^"' [Ln eK)<- inf I'^^"^' (//) . (3.42) 

(iii) Exponential tightness 

The sequence of distributions of Ln under P^''™' is exponentially tight. 

The proof of (iii) is given in Lemma 3.3 in Subsection 3.4. The proof of the lower bound (3.41) 

and upper bound (3.42) follows below. We are going to estimate the probabilities directly. The key 
idea is to split the two random mechanism. First, rewrite the sum over permutations in a sum over 
discrete empirical shape measures for integer partitions. The asymptotic of the discrete empirical 
shape measure is governed by the functional 5 (2.10). Any discrete empirical shape measure Q e Mn 
determines a probability measure P^' on (respectively on J7) , which represents the concatenations 
of Brownian motions with time horizon [0, (3] to Brownian bridges with time horizon an integer multiple 
of [0, P] depending on the length of the cycle present for the discrete empirical shape measure Q. There 
is a large deviations principle for Ln under Pj^^ with rate function 7^''' (2.12). 

Proof of the lower bound (3.41) of Theorem 2.2. Let G C M.i{0.) open with respect to 
the weak topology for probability measure. As $7 is Polish, the Levy metric is compatible with this 

topology ([DSOl]). Replace the sum over permutations in the definition (1.1) of the symmetrised 
measure P^^'"' by a sum over partitions in Vn- As in the proof of Theorem 2.1 in Subsection 3.1 we 
regroup for a given partition A G Vn the product of the Brownian bridges. In the following we write 
Tfe for rfc(A). Recall the definition of the measure fik = Mo,'o^ (3-26), and define for any partition 
A G Vn the measure (S^fcLi JPf"^* € M.i{^^), where Pfe is the normalised version of the measure jj-k, 
i.e., 

Ffe = /Xfe/Mfe = fik/Pk(3,A^ (0,0). (3.43) 

Hence, 

^N - z<^y^)^g. 2^ {II l^r, j^'^k ' 

N \t-'J Ae-Pjv k=l k=l 

where the volume factor in the numerator comes from the volume factor in the definition of the measure 



yLifc in (3.26). Thus, 
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where 

N 

j(g)ArQ(A;) 



pW)=(g)P^^«WGA^i(J^|^), QeMN. (3.44) 



k=l 



We now prove the lower bound for the large deviations principle in Theorem 2.2. o L^^^ is then 
a distribution in Miip.), i.e., we conceive Pj^' as measure on Q. 

Stirling's formula gives with the lower bound (3.31) the following lower bound 



^'■^"\Ln G G) > (1 - ^-dN/4p^N ^-^N"\ogN ^-Nlog{l+0{N)) 



^ Nd(k) ^ 



=(l_e-^W)iv^-iiv<^iogiv^-iviog(i+o(iV)) 1^ ^ exp(-iV5(Q))pj^>(L^GG), 

(3.45) 

where a G (^, 1), and where dist(0, 5Ajv) = ciV for some c > (see Proof of Theorem 2.1). 

We prove the large deviations lower bound for a sequence of boxes Ajv C M'^ with A^v T 1^*^ 
N ^ oo and dist(0, ^A^v) = cN for some c > 0. 

Proposition 3.2. Let (Ajv)iv>i be a sequence of centred cubes with Ajv T o-iT-d N/\A]\f\ p as 
N ^ oo and dist(0, d) = cN for some c > 0. Then 

liminf 1 log ( V exp ( - iV5(Q))p^H^iV G G) > - inf /(-"'(/z). 

QeMn 

Proof of Proposition 3.2. It suffices to construct, for any ^ £ G and Q £ Ai, some sequence 
iQN)N>i with Qpf G Aijsi and with ^ Q as N ^ oo strongly such that 

limsup (s{Qn) - ^ logP(«^)(Ljv G G)) < S{Q) + I^'^Hm)- (3-46) 

The first part of (3.46), i.e., limsupjv^go 5((5a) — '^iQ) is shown in the proof of Proposition 3.1 in 
(3.32). By d we denote the Levy metric on M.i{^), which generates the weak topology; see (3.47) 
below. 

By dist(;U, G) = infj^g^ d(/i, i') we denote the distance to a set G C Aii{C). Recall the Levy metric d 
on the Polish space A^i(O) [DSOl], defined for any two probability measures /x, G Mi{^) as 

d{n, u) = inf {(5 > 0: //(F) < u{r^) + S and u{r) < n{r^) + (5 for all F = f C O}, (3.47) 

where = {/j, e M.i{^): dist(/x, F) < 6} is the closed (5-neighbourhood of F. In the following let 
5iy such that Sjy I as N ^ oo. Let be a (5iv-neighbourhood of fj, in G, i.e. Ggj^ = {u E 
Mi{n): dist(zv,^) < 6n}- 

Let Qn & Mn defined as in (3.33) . We are going to use the Gartner-Ellis Theorem to deduce that 

liminf ^ log r«^)(-^^iv G G5^) > 

N—^00 N 
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For doing this, we first introduce for any F € Cb(i^) the logarithmic moment generating function 

TV 



Aiv(i^) :=logM^c«.)(e-(-'^^)) =log(nEo^fo(e-(-))"^^^'^) 
^ k=i 

N 

X^A.(/c)logEo^j(e^(^)), 



N 

n' 

k=l 



where E (q^) denotes the expectation with respect to the probabihty measure , and Eqq the 

expectation with respect to the probability measure P^. Here is supported on the path set (2.22). 
From the strong convergence of Qn to Q as ^ oo it is easily seen that A(F) := lim^v^oo jf^N{F) 
exists, and 

A(F) = ^Q(fe)logESj(e^(^)). 
fceN 

Since it is easily seen that A is lower semi continuous and Gateaux differentiable, and by the exponential 
tightness of the family P^^' o (see Lemma 3.4 in Subsection 3.4), [DZ98, 4.5.27] implies that 

liminf llogP^^'(L^ G GsJ > -I'^'iti). 

N—KX> iV 

This shows the second half of (3.46) and finishes the proof. □ 

Proof of the upper bound (3.42) of Theorem 2.2. From the exponential tightness (iii) in 
Lemma 3.4 it follows that there is a sequence of compact sets Kl C Ms{^) such that 



lim limsup^logf sup ¥'-^\L n e Kl)) = -oo. 



Hence, for the proof of the upper bound of the large deviations principle it suffices to show (3.42) for 
compact sets. Let K C Mi{^) compact. Then (3.30) and (3.44) imply 

P(7-'(L^ G X) < C7^"e^i°s(^+0W) J2 e^p[-NS{Q))F^^\Lj,eK) (343) 

QeMn 



We consider the logarithmic moment generating function of the distribution of Ln under P^\ 

N 

A^iF) :=logEp(a^)(e^<^'^-)) =iV^Q^(fc)logESj(e^(^)) ,FeCk{n). 



k=l 



Let now Qat G at be maximal for 

Q^e-^'5(^)p^^'"'(Ljv G K). 

The sequence {Qn)n>i has a limit Q G as ^ 00. Certainly, the limit A.{F) := limTv^oo jf^N{F) 
exists, and is lower semi continuous and Gateaux differentiable. Observe that I**^' is the Fenchel- 
Legendre transform of A. Now the Gartner-Ellis Theorem yields that 

limsup^P^^^(Ljv eK)<- inf I^'^\n). 

The cardinality of A^^v is bounded by (3.36), and by the continuity of Q 1-^ S(Q), the assertion follows. 

□ 
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3.4. Exponential tightness. In this subsection, wc prove the necessary exponential tightness asser- 
tions for the sequence of the empirical path measures under the symmetrised measures, and 
under the mixed product measures, Pj^\ for any Q ^ M.n- The proof of the latter exponential tight- 
ness is a variant of the standard proof for laws of empirical measures. Here, the main ingredient is the 
product structure of the probability measure. We equip the path space O with the following complete 
metric 

oo 

d(a;,e) = ^2-Mfc(a;,0, (3.49) 
k=l 

where 

Then d is a complete metric on fi, which generates the uniform convergence on compacts, and (O, d) 
is a Polish space. 

Lemma 3.3. The family of distributions of the empirical path measures Ljv under the symmetrised 
measure P^^™' is exponentially tight. 



Proof. The proof is in the spirit of the proof of [DZ98, Lemma 6.2.6]. Choose 5i,l G N, so small 
that for all A; G N, 

Pfof sup \Bs - Bt\ >\)< e-^\ (3.50) 

' ^ s,tg[0,fe/3], I ' 

\s-t\<ii 

Consider 



^/ = jo; G ri: u;(0) = 0,VA; G N: sup |a;s - a;t| < y|. 

\s-t\<h 



According to Arzela-Ascoli's Theorem, Ai is relative compact in O, because the set {ct;(0) : a; G A^} 
of the starting points is bounded, and for any A; G N we have 

lim sup sup \jJs — 0Jt\ = ^■ 

'-*oOajGA; s,te[0,fc/3], 

\s-t\<&i 

Put Ml := {n £ A4i(ri): /i(^) < j} and note that Mi is closed by Portmanteau's theorem. Let 
L G N be given and consider Kl H/^l ^l- ^asy to see that Kl is tight, hence Kl is compact 
by Prohorov's theorem. We shah show that F^^'^^Ln G kI) < e'^^ for any N eN, which implies 
the assertion. Observe that 

{Ln G M[} C {Hi G {1, . . . , A^}: (^p P^^ G > -}. 

Then, using the upper bounds in (3.30) and (3.37), we get 
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¥%''^\LMeMn<f'-r'U\ sup |i?«-5«i>y|>^) 

s-t|<5j 

-Ar5(Q) ^ 



/c{i,...,Af}, qeMn fc6{i.-,'v}, i';!ti<£ 

l^l>^v/( fcArQ(fc)>|/| ' 

< V V — TT e-' ^^^^^'^^ 

/c{i,...,Jv}, QeMn ^''^^ fce{i,...,JV}, 

l^l>^/' fcArQ(fc)>|7| 

< (^^2rAr«i(i-a) _ r^^^rAT"! ^ ^ ^^0(Ar)2Arg-Ar/_ 

JC{l,...,iV}, 
|J|>JV/( 

Hence, 

oo oo 
l=L l=L 

for all large N if R is chosen with 

(jj^OiN)2N+l < gArL(l-l/R)_ 

This ends the proof. 

□ 

Now we prove the exponential tightness of the empirical path measures under the probability 
measures P^',(5 G Mn, introduced in (3.44). 

Lemma 3.4. Let {Qn)n>i be a sequence of measures Qn G AIat such that Qn Q & M strongly 
as N ^ oo. Then the family of distributions of the empirical path measures Ljv under the measures 
Pj^^^ is exponentially tight. 



Proof. As in the proof of Lemma 3.3 there exists I G N and compact sets C such that for all 
keN 

Ff,{B ®kp G r^) < (e' - 1). (3.52) 

The set Ml = {u E i^(r^) < 1/1} is closed by Portmanteau's theorem. For L G N define 

= ni^L-^i- By Prohorov's theorem, each Kl is a relative compact subset of M.i{^). Then 
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Chebycheff's inequality gives that for any N G'N, and any Qn & Mn 



' JV 

N 



= e-^^'E^(^«,) (exp {2f 0, P^^' G r^})) 

fc=i 

fc=i 

N 

<e-2^' j;[(e')^'3'-(*^) <e-^'. 



fe=i 



Therefore, 



Pj^^'(L^ ^ i^L) < J]Pj^^'(L;v ^ Ml) < ^e-^' < 2e-^^, (3.54) 

l=L l=L 

which implies the exponential tightness. □ 

4. Appendix 

4.1. Large deviations. For the convenience of the reader, we repeat the notion of a large deviations 
principle and of the most important facts that are used in the present paper. See [DZ98] for a 
comprehensive treatment of this theory. 

Let X denote a topological vector space. A lower semi-continuous function I : X ^ [0, 00] is called a 
rate function if / is not identical 00 and has compact level sets, i.e., if /~^([0, c]) = {x ^ X : I{x) < c} 
is compact for any c > 0. A sequence (Xjv)iV6N of -Y-valued random variables Xn satisfies the 
large-deviation upper bound with speed ajv and rate function I if, for any closed subset F of X, 

hm sup — log ¥{Xn e F) < - inf I{x) , (4.55) 

and it satisfies the large- deviation lower bound if, for any open subset G of A:", 

liminf — logP(Xjv eG)<- inf I{x). (4.56) 

AT— >oo Ojv x&G 

If both, upper and lower bound, are satisfied, one says that (Xjv)jv satisfies a large-deviation principle. 
The principle is called weak if the upper bound in (4.55) holds only for compact sets F. A weak 
principle can be strengthened to a full one by showing that the sequence of distributions of is 
exponentially tight, i.e., if for any L > there is a compact subset Kl of X such that f{XN € KfJ < 
e-^^ for any N eN. 

All the above is usually stated for probability measures P only, but the notion easily extends to sub- 
probability measures P = Pjv depending on N. Indeed, first observe that the situation is not changed 
if P depends on A^, since a large deviation principle depends only on distributions. Furthermore, the 
connection between probability distributions P^v and sub-probability measures Pjv is provided by the 
transformed measure Piv(^Ar & A) = ¥]\f{X]\f G A)/¥n{Xn £ X): if the measures Pjv o X]^^ satisfy a 
large deviation principle with rate function /, then the probability measures P^v o X^^ satisfy a large 
deviation principle with rate function / — inf /. 

One standard situation in which a large deviation principle holds is the case where P is a proba- 
bility measure, and Xn = ( Yi -|- • • • -|- Y^) is the mean of N i.i.d. -Y-valued random variables Yi 
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whose moment generating function M(F) = J e^^"^^^ dP is finite for all elements F of the topologi- 
cal dual space X* of X. In this case, the abstract Cramer theorem provides a weak large deviation 
principle for {X]\f)N^fq with rate function equal to the Legendre-Fenchel transform of logM, i.e., 
I{x) = su-pp^;t;*{F{x) — logM(F)). Au extension to independent, but not necessarily identically 
distributed random variables is provided by the abstract Gartner-Ellis theorem. 

In our large deviations results we shall rely on the following conventions. For X = Q or X = 
fifc. A; G N, we conceive Aii{X) as a closed convex subset of the space X = Ai{X) of all finite signed 
Borel measures on X. This is a topological HausdorfF vector space whose topology is induced by 
the set Cb(^) of all continuous bounded functions X ^ M. Then C]3{X) is the topological dual of 
Ai{X) [DSOl, Lemma 3.2.3]. When we speak of a large deviation principle for A^i(X)-valued random 
variables, then we mean a principle on Ai{X) with a rate function that is tacitly extended from 
Mi{X) to M{X) with the value +00. 



4.2. Bose functions. These functions are defined by 



9s{a) 



r{s) 



— dt = y k~'e-'^'' for all a > and all s > 0, 



(4.57) 



k=l 



and also a = and s > 1. In the latter case. 



k=l 



which is the zeta function of Riemann. The behaviour of the Bose functions about a = is given by 



9s{oi) 



r(i-,)a-i + Er=oC(5-^) 



i-ar 



{s-iy. 



k\ ' 



s^l,2,3,... 



k^s-l 



At a = 0, Qsia) diverges for s < 1; indeed for all s there is some kind of singularity at a = 0, such as 
a branch point. For further details see [Gra25]. 



4.3. Quantum statistics for Bosons. We review some basic facts on quantum systems of non- 
interacting Bosons. Fix a box A C W^. We consider a system of A'^ identical Bosons enclosed in A 
in thermal equilibrium at inverse temperature (3 > 0. The total energy of the system is given by a 
Hamilton operator 

N 

1=1 

where — Aj is the Laplacian with some boundary conditions representing the kinetic energy of particle 
i. The states of the system belong to the symmetrised subspace H^^'''^'^^ of the A/"-fold tensor product 
nf^ of the one particle Hilbert space Ha = L^(A). The projection P^'""* oiUf^ onto U'-^''"'"^^ reads 

p(sym) ^ \ ^ TT 

where : is a unitary representation of the permutation group 6;v on defined by 

C4(V'i <^ • • • <^V'Ar) = V'Tra) • • • ® ^7r(Ar), V'fe ^ T^A, k = 1, . . . , N,Tr E &N- 
The quantum canonical partition function is 

ZAiP,N) = TV^^;v,...)(e-^<') = TV^^.(pj^-)e-/^<') = ^ E /aW, 
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where /a is a function on © at defined by 

/A(7r) = Tr^^N (c/^e-^^i^') , Tre&N- 

Each permutation tt G ©at can be written as a product of independent cycles. A cycle of length k with 
1 < < A?" is a chain of permutations, such as 1 goes to 2, 2 goes to 3, etc. until A; — 1 goes to k and 
finally k goes to 1. A permutation vr G Sjv may be built up of ri 1-cycles, r2 2-cycles,. . . , ttv A-cycles; 
we then say tt has the cycle structure or is o type (ri, . . . ,riv) with the X^fcLi ^i^k = A"- Hence, any 
partition of the integer A" defines a conjugacy class of permutations such that the trace operation for 
the elements of this class remains constant. The genuine task of quantum statistical mechanics is to 
prove and analyse the thermodynamic limit lim^_^]jd 1/|A| log Z\{p, N) such that A/|A| p £ (0, oo) 
([Rue69]). Feynman 1953 introduced the functional integration methods for traces, see [Gin 70] for 
details. Calculations in the Fourierspace and in the grandcanonical ensemble go back to Bose and 
Einstein 1925. Bose-Einstein condensation is defined as a macroscopic occupation of the zero mode (in 
Fourierspace), which corresponds to particles with zero wave vector, i.e, zero momentum. This is only 
present for systems of Bosons with given density at sufficiently low temperatures. A mathematical 
criterion for Bose-Einstein condensation has been proposed by Onsager and Penrose [OP56]. Bose- 
Einstein condensation for non-interacting particles is equivalent to the appearance of long cycles 
(cycles that grow as the system size) ([Sii02],[Uel06]). Proofs for Bose-Einstein condensation in the 
thermodynamic limit have been obtained for mean-field models ([BCMP05],[DMP05]), and for lattice 
systems with hard-core exclusion interaction in [LSSY05]. For systems confined in trap potentials, 
Bose-Einstein condensation was proved in the so-called Gross-Pitaevskii limit, which is different to 
the usual thermodynamic limit ([LSSY05]). 

Acknowledgments. The authors thanks Tony Dorlas and Joe Pule for discussions on cycles at 
Dublin Institute for Advanced Studies. 
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